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Abstract 

The temperature renormalization group equation (TRGE) is compared with a diagram- 
matic expansion for the (0 4 )4-theory. It is found that the one- loop TRGE resums the 
leading powers of temperature for the effective mass. A two-loop contribution to TRGE 
is required to do the leading resummation for the coupling constant. It is also shown that 
the higher order TRGE resums subleading powers of temperature. 
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1 Introduction 



The temperature renormalization group equation [|TJ (TRGE) has been used by several 
authors to study the high temperature behaviour of different theories [H, fH, H, |5|, H • With 
few exceptions it has not been clear exactly what kind of improvement the TRGE 
implies. In some cases [fj, || a 1-loop approximation of the effective potential was used to 
determine the renormalization flow of the mass and coupling constant at zero momentum. 
This gives the first correction in h but not necessarily the leading correction in T/m as 
we shall see. The purpose of this paper is to give a more detailed discussion of the TRGE 
for the (0 4 )4-theory. 

When calculating thermal corrections perturbatively in a QFT at finite temperature it 
is often found that the 1-loop result goes like a power of T at high temperature. The 
correction is, therefore, always large for large enough T which makes the perturbation ap- 
proach unreliable. The self-energy, for instance, is the correction to the effective (mass) 2 
of a particle surrounded by a heat-bath. It gets a correction going like T 2 at high T 
in the (0 4 )4-theory. If the quasi-particles in the thermal environment have a mass very 
different from the zero temperature one, it is clear that it is not a good approximation to 
expand around the zero temperature excitations. On the other hand, expanding around 
quasi-particles at one temperature to calculate the self-energy at a temperature that is 
close has a better chance to succeed. Also, a sequence of many small steps in temperature 
is intuitively more likely to give a better result than taking only one big step. The situ- 
ation is familiar from zero temperature calculations where the mass can be renormalized 
at different scales fi. Using a mass defined at fi to calculate the mass at the scale // gives 
rise to log(/x'/ /z) terms which can be large. There the remedy is well-known as the renor- 
malization group which sums up small steps in \x. The renormalization group equation 
(RGE) at T = can be formally derived from the freedom of choosing fi in the renormal- 
ization prescription. At finite temperature, the mass of the quasi-particles, around which 
to expand, can be defined at any temperature and there is a freedom of choosing that 
temperature which should not influence on the final answer. This freedom can be used to 
derive the TRGE which implements the summation over small temperature steps |l]. 

The result from solving the 1-loop TRGE is certainly different from just computing a 
1-loop correction but it is not obvious exactly what the difference is. It is known that 
a resummation of leading log/z terms is obtained by solving the RGE at zero temper- 
ature. Since there is a similarity between log/i and T it is natural to guess that the 
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TRGE resums the leading powers of T. This is true for the (0 4 )4-theory and that is the 
basic result of this paper. It turns out, however, that the TRGE can only resum certain 
subdiagrams in an effective way. By effective I mean getting a result to the n-th order 
without doing an n-th order calculation. The diagrams that are effectively resummed are 
the ones containing loops that give a leading power of T 2 (hard thermal loops). These 
terms can also be resummed using an effective propagator with a T dependent mass f7|. 



The relation between the two methods is discussed in Section 3.1 



In Section ^| the thermal 9 and (3 functions are discussed and their relation to the 
diagrammatic expansion of m(T) and g(T). The TRGE to leading order is solved in 
Section WA. 



2 Resummation of leading powers of T 
2.1 General ideas 

Our goal is to calculate the mass and the coupling as functions of the temperature 
(m(T), g(T)). The meaning of m and g depends on the renormalization conditions (RC). 
To give a physical interpretation of m(T) we define it to be the real part of the pole of 
the propagator of physical excitations. This assumes that it is meaningful to talk about 
single stable excitations which may be difficult at finite temperature when the excitations 
eventually dissipate away, but we take it as an approximation. Adding an imaginary part 
to the mass would take into account the dissipation. 2 In some cases it may be needed to 
to consider excitations with a spectrum that differs drastically from the free one, but we 



do not discuss that here (see Section |3^2|) . The coupling constant is related to experiment 
through some scattering process and we define it to be the real part of the 4-point function 
with external momenta corresponding to that scattering. For simplicity we assume that 
the momentum depends only on the mass and some parameters that are kept fix, such as 
scattering angles. The renormalization point is then denoted by pm = Pi{w,R)- Thus the 
renormalization condition is 

Rer( 2 )(p,T)L 2= , t=Tr = 0, 



2 Perturbative propagators with complex poles do not admit any spectral representation, which the 
full propagator can be shown non-perturbatively to have (at least at equilibrium). It may therefore be 
necessary to consider some generalization of the TRGE to the spectral function if dissipation is to be 
taken into account consistently. I want to thank Dr. Umczawa for pointing out this to me. 
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ReT {4) (pi,T)\ Pi{mR)jT=TR = -g R (T R ) . (1) 



jReT^\ p , =m 2 wT=TR = -l , (2) 



In addition, the wave function has to be renormalized and we choose 

_d 

dp 2 

as RC. Other conditions (like the MS scheme) can be imposed if the only purpose is to 
eliminate the UV divergences. Here we also want to give the expansion parameters m R 
and g R a physical meaning. 

From the RC the corresponding 9, ft and 7 functions can be defined ||TJ 

1 d 

Ri m Ri Tr) = ~^^Qf ReTi2) (P^ T ">\ T R ' 



ft(g R ,m R ,T R ) = ^g R --^ReV^(p iR ,T)\ TR 



1 d d 

l{g R ,m R ,T R ) = -2 9T^ Rer(2) ^' T ^ Pfl ' Tr ' ^ 
The mass and coupling, defined through Eq.([TJ), are obtained by integrating 9 and ft, 

% = 0(g,m,T) => \ m = m(T;g R ,m R ,T R ) 
% = ft(g,m,T) \ 9 = g{T;g R ,m R ,T R ) 



and we get in this way the zero of Rer^(p, T) and the value of ReT^ (p R , T) as functions 
of g R , m R , Tr and T. 

The same functions could in principle be calculated by the diagrammatic perturbation 
series as an expansion in g R . The idea here is to compute 9 and ft perturbatively and 
then solve the non-linear ordinary differential equations Eq. @ which are non-perturbative 
relations. To study the relation between the two approaches it is convenient to use the 
invariance under a change of the initial conditions to write the linear partial differential 
equations 

( d d d \ 

[t^t + 0(g R ,m R ,T R )- \- ft(g R , m R , T R )— — m(T; g R ,m R ,T R ) = 0, 

\ a R om R og R 1 



d d d \ 

+ 0(g R ,m R ,T R )- h ft{g R ,m R ,T R )—— g(T; g R ,m R ,T R ) = 0. (5) 



oT R dm R dg R/ 



The definition of 8, (3 and 7 differs from Ref. M by factors of Tr and mR. 
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I(T,m R ) - I(T R ,m R ) 
Figure 1: The tadpole plus its counterterm. 



If now 9 and (3 are computed to finite order in g R what do the corresponding solutions 
for m and g mean in terms of perturbation theory? To answer that question we make 
an Ansatz for the solution as a power series in g R and determine the coefficients. Before 
doing that we shall take a closer look at the leading T behaviour of 9 and (3. We must 
also be aware that, by making the Ansatz above, we assume a power series expansion of 
m and g around g R = 0. 

2.2 Power series for 9 and (3 

We are interested in the leading T dependence of 9 and (3 at each order in g R , so we have 
to find the diagrams that contribute most. According to Eq.(|3|) the 9 and (3 functions 
are obtained by taking suitable T and p derivatives of Rel^ 2 ) and Rel^ 4 ) and putting 
p = p R , T = T r afterwards. The Rer 1 - 2 ) and Rer*- 4 - 1 should be renormalized with the 
physical self-consistent RC in Eqs. (|l], 0) . This turns to be crucial. 

The dominant diagram at high T is the tadpole correction to the mass which gives 
a factor T 2 . All other loops give a factor T. Therefore, the dominant diagrams to each 
order in g R are the ones with as many tadpole subdiagrams as possible. This is true for 
However, for each such tadpole there is a corresponding diagram with a counterterm 
defined at T R (see Fig.QI|)). 

Let us denote the tadpole contribution by I(T,m R ). The sum of all diagrams with 
n tadpoles and their counterterms contains a factor (I(T,m R ) — I(T R ,m R )) n . Evidently, 
after taking the the T derivative and putting T = T R , this factor vanishes for all n > 2. 
We conclude that only diagrams with at most one tadpole as subdiagram contributes to 9 
and (3 and thus the n-th order term in for 9 and (3 goes like g n R T n and g R +1 T n respectively. 
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We write the high T expansion as 

rri OO rji OO rp 

P(— ) = ^E4(— ) fc E(— 




E^) fc E(-^r p ^- 

k=i m R P =o m R 



(6) 



where k is the order to which 9 and /3 are calculated. By looking at low order diagrams 
we see that /3io is zero which is important in next section. 

In addition to the terms in Eq.(^) there are subleading terms like T m ln n T. If they 
are considered to be of the same order as T m independently of n, they do not affect the 
ordering of terms with different powers of T that is used in next section. They also do 
not occur to the leading order. Therefore, I do not write them out since they would only 
complicate the formulas. 

2.3 Power series for m and g 

The purpose of this section is to relate the solutions of Eq.(||) to the perturbative calcu- 
lation of 9 and (3. A diagrammatic expansion of m and g yields a power series in gji so 
we make the Ansatz 



why we related pn to m,R. it reduces the number of dimensionful variables. Of course, 
we could have related pn to Tr and get the same simplification here but with an other 
physical interpretation of m and g. The leading behaviour of 9 and j3 in last section 
was considered for a pr which does not depend on T R . We could also have treated the 
renormalization group equation in both pr and Tr as proposed in . 

Our primary interest is to see how the T dependence in 9 and (3 gives a T dependence in 
m and g so we further expand M n and G n 





T T, 



R 



E( — ) g ( — ) l Mn ql 

^ m R m R 



m R rriR 



T T, 



R 




(8) 



nsu i 



m R m R 
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where the summation range for q, I, s and u is yet to be determined. The comment at 
the end of last section regarding possible In T factors apply to this expansion as well. 
From Eq.(||) and Eq.(^) we find 

n oo 

qM nql = M n _ k ^ 1+p _ k j{6 kp (2 + k-p-q-l) + f3 kp (n - k)] , 

k=l P=o 



sG nsu = -Y,Y. G n-k,s-i+p-k,u[^kp{^ + k-p~s-u)+f3 kp {n-k + l)} , (9) 

k =lp=o 

by identifying the coefficients of the same powers of g R , T R and T. The Eq.@ determines 
M nq i (G nsu ) in terms of M n i q i (G n / su ) where n' < n, and it can be applied inductively 
except for q = (s = 0) 0. 

We are interested in the T dependence but Eq.(H) is a differential equation in T R . The 
T dependence enters through the boundary conditions which also determines M nqi {G nsu ) 
for q = (s = 0). The boundary conditions are 

m(T R ;g R ,m R ,T R ) = m R , 
g(T R ;g R ,m R ,T R ) = g R , 

or in terms of M n and G n 

M n {^-- ^) = , n > 1 ; M = 1 , 
m R m R 



G n (— ;— ) = 0, n>l; G = 1 . (11) 



They lead to 

M n0 Z = -ErM^j^ , 

where X)' means that r = (v = 0) is excluded. 

We shall now determine the summation ranges for q, I and s,u in Eq.([8]) using Eq.(Q) 
and Eq . flT2|) . It is obvious that there is no lower limit in general. We define the upper 
limit by 

Q(n, I) = max q for which M nql ^ 
L(n, q) = max I for which M nq i ^ , 



Again, when including InT terms this exception is modified. 
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(14) 



S(n, u) = max s for which G nsu ^ 

U(n, s) = max u for which G nsu ^ . 
When doing this we must be careful with the last factors in brackets in Eq.(^) which can 
be equal to zero for certain values of k,p,n,q,l,s and u (remember also that p w = 0). 
As noted before, their form is different when including possible InT factors. The easiest 
way to find Q, L, S and U is to do an explicit calculation for small n, guess the general 
form and then prove it by induction. There are some comments about the procedure in 
Appendix |A|. 



The initial condition for n = is 

Q(0,0) = 



£(o,o) 

5(0,0) 
*7(0,0) 









Q(0,Z ^ 0) = -oo 
L(0,q ^ 0) = -oo 
S (0, it 7^ 0) = -oo 
U(0,s^ 0) = -oo 



The result for n > 1 is 



Q(n,l<0) = 2n-2 
Q(n,0 <l <2n) = 2(n - [ 
Q(n, I > 2n) = — oo , 

L(n,q<0) = 2n-2 
L(n,0 < q <2n) = 2(n 
L(n, q > 2n) = — oo , 

S{n,u<0) = 2n-2 

S(n,u = 0) = 2n-l 

S(n, 1 < u < 2n - 1) = 2(n - 1 

S(n, u > 2n — 1) = — oo , 

U(n, s < 0) = 2n - 3 

[7(n, s = 0) = 2n - 1 

C/(n, 1 < s < 2n - 1) = 2(n - 

t/(n, s > 2n — 1) = — oo 



¥1 



s+l 
2 



]) 



(15) 



(16) 



(17) 



18) 



(19) 



Here, [q] means the largest integer less or equal to q. These upper limits show that the 
leading T dependence for each order in g R goes like 



/;/ 



oc m R J2g n R T 2n 



n=0 

oo 

g oc toE^T 2 "" 1 

n=0 



(20) 
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which could be seen directly from the diagrammatic expansion. We shall now find out 
what values of k that enters in the determination of the z-th subleading power of T in M n 
and G n . Therefore, we look at the terms in Eq.(^J) with I = 2n — i and u = 2n — 1 — i, 
which are given by 

E(~)' ?M «,9:2n-i = - E ( — Y E E -M n - Kq - 1+p - K2 n-i 

q " LR q<Q(n,2n-i) " lR fc=l P=0 ^ 

2n—i<L(n,q) 

[6 kp (2 + k-p-q-2n + i) + (3 kp (n - k)] 



n oo ^ 

+ X ' EE - M n-k,r-i+ P -k,2n-i-r[0k P (2 + k-p-2n-i) + (3 kp (n - k)] 

r<Q(n,2n-i-r) fc = l p=0 ^ 

2n— i<L(n,r)+r 



(21) 



rp rp n OO 

E(~) SG «.^2n-l-i = - E ' (~) S E E -G n -k,s-l+p-k,2n-l-i 

s m R .<s(«,a»-i-*) m « fc=lp=0 S 

2n-l-j<C/(n,s) 



[# fcp (2 + A;-p-s-2n + i)+ /3 fep (n - fc + 1)] 



n oo 

+ E ' EE -G n -fe ir _i +p _fc,2n-i-i-r[#fcp(2 + A;-p-2n-i)+ (3 kp (n -k + 1)} 

r<S(n,2n-l-i-r) fc = l p=0 ^ 

2n-l-i<U(n,r)+r 



(22) 

The first term on the RHS is for q ^ (s ^ 0) and the second is for q = (s = 0). The 
general solution to the problem of finding what values of k contribute for a given value of 
i is rather complicated because of the many special cases that have to be considered. The 
cases of i — 0, 1 are treated in next section. The way of solving the problem for general 
i is discussed in Appendix B, here I just state the result. To resum the i-th subleading 
powers of T using the TRGE, the 6 and (3 functions have to be calculated to the fc-th 
order where k satisfies 

m : k < i/2, all p , . 

~ 1 ' y (23) 
oip + k<i + l 
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g : k < i/2, all p 
orp + k<i + 2 



(24) 



Remember that the % is defined as T 2n ~ % for m and as T 2 "^ 1 ^ 1 for (7. 

A discussion of the result is given in Section ||] and here we just notice that the solution 
of the TRGE, using a perturbatively calculated 9 and (3 (i.e. finite k), indeed results in a 
resummation of leading and sub leading factors of — . 

ran 

2.4 Special cases when i = 0, 1 

The most important cases are, of course, the leading and first subleading contributions. 
We start with the leading terms for the mass which go like in Eq.(p0|). The analysis in 
Appendix [B] gives that only k = 1, p = contribute, i.e. 8\o- That can also be seen from 
a diagrammatic expansion of m 2 ; only the 1-loop tadpole give a gRT 2 dependence, all 
other diagrams give at most g R T 2n ~ x . Therefore, the leading T dependence of m(T) is 
given by 

m(T) = ^m R + g R 6 w {T 2 -T 2 R ) . (25) 

The next-to-leading order (z = 1) requires the knowledge of #10, 6*n, #20, Pu and /?2o (from 
p + k < 2, Eq.([23|)) which gives a much more complicated TRGE to solve. 

Let us continue with the expansion of g(T) which gives a more interesting result. It 

turns out that to leading order (i = 0) it is necessary to calculate 

0io, #u and P20 
5 

That is, to get the leading T dependence we must actually do a 2-loop calculation! 
Since the leading order TRGE can be solved explicitly we shall take a closer look at it to 
see what the actual values of 9, (3, m{T) and g(T) are. What we need to know is #10, fin 
and /?2o, and they are obtained from the diagrams in Fig.(Q), where we use the notation 

r°° dk k^ 1 

F^T,m)=l — r -^ JT — i ,u; = Vk 2 -T^. (26) 
From these diagrams we find the TRGE 

dm gT 

= 24^' { ' 

5 From Eq.(p4|) we havep+k < 2, but On and 620 do not contribute since the bracket [2+k— p— s~ 2n+i] 
vanishes, see Eq. (]22]) . If we had included the InT terms the bracket would in general not be of the same 
form. However, to highest order there are no InT terms so the conclusion is correct. 
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9R T?2lrn „ \ 9RT 2 



;F?(T,m R ) 




4vr 2 1V ' n/ 24 
3g R dF?(T,m R ) _ 3g R T 




An 2 dm R lQnm R 




16tt 4 dm% 2Xitm% 



3,J W.m^ 9rTTx 



16tt 4 dmfj 1V 2567rm^ 

Figure 2: Diagrams contributing to the lowest order 6 and f3. The vertex is renormalized 
at zero momentum. 

= ^ I 9 " T2 (28) 

dT IQnm 1287rm 3 ' 1 ; 

The last equation can be written as an exact differential and the solution when T R = is 

9(T) = 1 + ^ • (29) 

16tt m(T) 

Then the equation for m(T) reads 

dm = 2ng R T 

dT 48vrm + 9^T ' 1 ' 

which, after a change of variables m(T) = t> (T) T, turns into a separable equations and 
can be solved by integration. The implicit solution is 

167r m 24 mR [g 67rm _ T( v /81^ + 384tt 2 (? r - 9^)/ 
In the T — > 00 limit we find that 



m(T) /gfl 3g fi 3 to 

~T~ V24 + W 32^' [6I) 
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and 



g(T) -> g R 



8lg R + 384^g R -9g R 



(33) 



/81^ + 3847r 2 ^ + 9^ 

In order to determine the subleading behaviour of g(T) we would need to calculate 6*10, 
020, P12, P20, P21 and fyo. Note that the subleading terms can become dominant 
since the leading terms sum up to a constant when T — ► 00 (see Section EO) . 



3 Conclusion and outlook 

The main conclusion of this paper is that the TRGE resums the leading and subleading 
powers of T in much the same way as the T = RGE resums leading In /1 terms. It is 
also interesting to notice that a 2-loop calculation is required to get the leading behaviour 
for the coupling constant. This was observed by Funakubo and Sakamoto in Ref. [[3] who 
studied the O(N) model in the limit N —>■ 00 where it is exactly solvable. 

It is, of course, not a miracle that some higher order diagrams are determined in terms 
of lower order diagrams. At zero temperature, the leading In \x contribution comes from 
diagrams with many divergent subdiagrams. It is those subdiagrams that are controlled 
by the f3 function. Higher order diagrams with leading ln/z terms are constructed by 
multiple insertion of those subdiagrams. At finite temperature there is one type of loop 
that differs from others: the quadratically divergent tadpole which gives a factor T 2 [7|, || . 
All other loops give a factor T (in the high T limit). At the n-th order, if we want the 
n-th subleading contribution to m(T) (oc g R T n ) all diagrams to that order contribute. 
It is, therefore, natural that we have to calculate the 9 and (3 functions to n-th order 
(actually, even if we want the (n — l)-th subleading term we get a k = n contribution 
from the condition p + k < i + 1 (Eq.(p3|)) when p = and i — n — 1). We can briefly say 
that the TRGE is a way of resumming the subdiagrams with a T 2 dependence of which 
there is only the 1-loop tadpole in the (0 4 )4-theory. 

A difference between the zero and finite T RGE is that the f3 function at zero T 
does not depend on \i in the MS scheme, while the (3 function at finite T depends on T. 
That makes it harder to prove a leading T resummation and it works only because the T 
dependence is not too strong, only g R T n and not g R T 2n . The reason for that is, as we saw 
in Section |2.2| , that the highest power of T is subtracted off by the counterterms. Then 
we are back to the physical starting point, namely to use the physical T dependent mass 
already in the propagator. In that way the highest T dependence does not contribute to 
the P function. 
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3.1 Effective propagators 

The method of using an effective propagator with a mass term m 2 R + g#T 2 /24 also resums 
the leading powers of T to each order in g R . By definition of the effective mass we have 
the same expression for m(T) when Tr = (Eq.(E5j)) and from the 1-loop diagram in 
Fig.(0) we get for the coupling constant 



g(T) = g R + ^- ln(l + ^fV) - ^ , Z (34) 



which agrees with Eq.(|29|) to leading order. 

Still they predict quantitatively different behaviours at very high temperature. The 
logarithmic term in Eq. (|34"D , that comes from the renormalization at T = 0, is formally of 
lower order in T/tur and should be neglected. Then Eq.(|3"4"D gives a negative effective 4- 
point coupling at sufficiently high T if qr > 327r 2 /27, which indicates an instability even 
though it can be stabilized by higher n-point functions (such as the 6-point function). 
On the other, hand the logarithmic term is actually dominant at high T and gives an 
increasing g(T). The TRGE result in Eq . (|33f) does not show any of these behaviours. The 
difference between the effective propagator approach and the TRGE is in the subleading 
terms, and it is clear that if the leading terms sum up to something finite when T — > oo 
the subleading terms can become dominant. That is also the reason why m(T) and g(T) 
are found to diverge at a critical temperature in Ref. 0] but not here. 

A problem with the effective propagator is that it leads to T dependent infinities in 
the tadpole diagram in Fig.(|I]) when the mass is explicitly T dependent. It occurs in 
subleading terms and is subtracted away if all diagrams to that order are included ||. 
Also, in this was pointed out as a defect of not treating subleading terms consistently 
while in |10| it was considered to be legitimate with such T dependent infinities. In my 



point of view the bare parameters of the theory has nothing to do with the thermal state 
of the system and must, therefore, be temperature independent. The TRGE is, however, 
derived from the T independence of the bare parameters and satisfy this criterion by 
definition. It is a subject for further study to understand significance of subleading terms 
in the TRGE approach. 



3.2 Other theories 

As we now understand, the method of resummation using the TRGE works because the 
leading temperature corrections can be accurately approximated by a T dependent mass 
term. In the (0 4 )4-theory it is particularly simple since the tadpole does not depend on 
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the external momentum. The possibility of extending the result of this paper to other 
theories is related to the possibility of including the leading T dependence in the propa- 
gator, i.e. to expand about the physical excitations at finite temperature. For instance, 
the (0 3 )6-theory has a momentum dependent 1-loop correction to the mass. It turns out, 
however, that the leading T 2 part of that correction is momentum independent and the 
TRGE should work even for that theory. 



In QCD, the gluon polarization tensor has a T 2 term which is momentum dependent 
and cannot be approximated by a constant 11, 12 1. This indicates that the TRGE, as 
described here, does not yield a consistent resummation. It may be possible to generalize 
the TRGE to include a momentum dependent mass term. The choice of RC determines 
the propagator and we have only considered finite changes in a momentum independent 
mass by renormalizing at a fixed momentum. In principle, one can allow for a momentum 
dependent mass term. The reason for not doing so is that it makes calculations more 
complicated, but if the physical excitations cannot be well approximated by free particles 
with T dependent mass, one has to use something better. A perturbation method using 
effective propagators and vertices to include the T 2 contributions was developed in j!3 



The analysis in this paper indicates that a 2-loop calculation would be necessary also in 
QCD to get a TRGE which resums the leading powers of T for the coupling constant. 



I want to thank Drs. Umezawa and Sakamoto for discussions and The Theoretical Physics 
Institute at University of Alberta for their hospitality during my stay. I also thank The 
Swedish Institute for financial support. 



A Appendix 

We out-line the determination of the functions Q,L,S, and U defined in Section |2~3j From 
Eq.(||) we see that 

Q(n, I) = max[Q(n - k, I) + 1 - p + k] (35) 

k,p 

where 1 < k < n and p > 0. If the result is less or equal to zero we must also check 
Eq . (jjjD which gives a lower limit 

Q{n,l) > if l< L(n,r) +r ; r ^ . (36) 
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Since the determination of Q(n,l) only involves Q(n',l) with n' < n we can solve the 
problem iteratively in n with initial condition 

Q(0,0) =0 , Q(0,/^0) = -oo . (37) 

Each time we determine k,p in Eq.(^) we must be sure that the factor in brackets in 
Eq.(||) does not vanish for these values. It does not pose any problem here but it does 
when calculating S and U . 

Similarly L is then determined by 

L(n,q^0) = max k , p [L(n- k,q- 1 +p- k)] ^ 
L(n, 0) = max r7 L [L(n, r) + r] , 

with the same initial condition as for Q. 

The procedure is the same for S and U but we find that 

5(1, 0) = max[£(0, 0) + l- p + l]=2, (39) 

occurs for k = 1, p = 0. This term is, however, absent in Eq.(|9|) since /3xo = an d 
1 + k—p — s — u = 0. Therefore, S^ljO) = 1 for k — 1, p — 1. That is the major difference 
between Q, L and S 1 , [/. As we can see from Eq.(§) S(n, 0) and U(n, 0) has to be treated 
separately from S(n, 1 < u < 2n — 1) and U(n, 1 < s < 2n — 1). 



B Appendix 



In order to determine which values of k and p that occurs in the summations in Eqs. (|2T| , 22[) 
we look at the first and second term separately. 



In the first term in Eq.fl21"|) the summation over q runs through all values satisfying 

q < Q(n,2n-i) 
2n — i < L(n, q) . 

For each such q the terms that contribute to the k and p summations satisfy 

q—l+p — k < Q(n — k,2n — i) 

2n — i < L(n — k, g — 1 + p — k) . 

Finally, the factors in brackets in Eqs.(|2l|,|22|) must be non-zero. This analysis gives that 
for i > 2n there is no restriction on k and p (other than 1 < k < n and p > 0, of course) 
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and for 2 < % < 2n we have k < i/2. There is no contribution from the first term when 
i = 0, 1. 



The summation range for the second term is determined by 

r < Q(n, 2n — i — r) 
2n — i < L(n, r) + r , 

and the terms in the k,p sums satisfy 

r—l+p—k < Q(n — k,2n — i — r) 
2n — i — r < Lin — k,r — 1 + p — k) 

We find that this term contributes whenever p + k < i + 1. 



(42) 



(43) 



In a similar way we check the summation in Eq.(|2"2"|) where in the first term the sum- 
mation runs over all s satisfying 

s < S(n,2n-l-i) 
2n-\-i < U(n, s) , 

and the k,p summation ranges over 

s — 1 + p — k < S(n — k, 2n — 1 — i) , , 

(45) 

2n — 1 — i < U(n — k, s — 1 + p — k) . 



We find, using Eqs.([T8|, [T9|), almost the same conditions as for M, i.e. k < i/2 and 
p + k < % + 1. However, the expression for S(0, u) and U(0, s) do not follow the gen- 
eral form and give exceptions for n = k in the RHS of Eq.(^). This gives an extra term 
when k = n and p < i + 2 — n. The highest n (=k) for which this term can contribute is 
when p = 0, i.e. n = k = i + 2. 



In summary, we have found that for each iteration of the solution of Eqs.(|2~T|,|2"2"D (if 
we imagine an iterative solution in n) the terms contributing are the ones satisfying 

k < i/2, all p (also 1 < k < n) , 

ItV--Ll " (46) 

or p + k < % + 1 . 



In addition, in the [i + 2)-th order in Eq . fl22|) there is an extra term with 

p + k < i + 2 . (47) 
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We conclude that, in order to get the i-th subleading term of the T^/m^-expansion of 
7n (T) we need to compute the 9 and (3 functions to the (i + l)-th order (respectively 
the (i + 2)-th order for g(T)). Actually, we only need the leading T-dependent part of 
the highest order diagram when calculating 9 and (3. The Eqs.(|46|,|47|) tell more precisely 
what is needed. An explicit example for i = is given in Section |2.4j . 
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